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ABSTRACT: We compute the — > e conversion in the type-I seesaw model, as a function 
of the right-handed neutrino mixings and masses. The results are compared with previous 
computations in the literature. We determine the definite predictions resulting for the ratios 
between the jx — > e conversion rate for a given nucleus and the rate of two other processes 
which also involve a fi — e flavour transition: \x — > e"f and fi —> eee. For a quasi-degenerate 
mass spectrum of right-handed neutrino masses -which is the most natural scenario leading 
to observable rates- those ratios depend only on the seesaw mass scale, offering a quite in- 
teresting testing ground. In the case of sterile neutrinos heavier than the electroweak scale, 
these ratios vanish typically for a mass scale of order a few TeV. Furthermore, the analysis 
performed here is also valid down to very light masses. It turns out that planned fj, — >• e 
conversion experiments would be sensitive to masses as low as 2 MeV. Taking into account 
other experimental constraints, we show that future fj, — > e conversion experiments will be 
fully relevant to detect or constrain sterile neutrino scenarios in the 2 GeV— 1000 TeV mass 
range. 
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1 Introduction 

The recent experimental evidence for neutrino masses has shown that lepton flavour is violated 
in the neutrino sector. This inevitably leads, at the tree or one-loop level, to rare processes 
violating charged lepton flavour, such as I — > I'^y, I — > 31' or /itoe conversion in atomic nuclei. 
Current experimental bounds on the various rates are expected to be improved in the near 
future by a long series of new experiments. In particular, fi — > e conversion processes [1] 
will become especially competitive, as the sensitivities for various nuclei are expected to be 
improved by several orders of magnitude, 
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as compared to the present sensitivities 



R^Ue < 4-3 x 10 
Rt» e < 7 x 10- ] 
R™ < 4-6 x 10 
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[9]. 



The rates predicted by neutrino mass models are in general expected to be far below 
these sensitivities. However this is not necessarily the case. The rates are highly model 
dependent. Classes of models, such as those based on approximate conservation of lepton 
number (predicting 2 or more quasi-degenerate heavy fermions), can naturally give large 
flavour-changing rates and the experimental sensitivities can be saturated without fine-tuning. 

The fj, — >■ e conversion rate R^ e has been calculated in the literature for the various 
possible types of seesaw, with right-handed neutrinos [10-16], scalar triplet(s) [17-19] and 
fermion triplets [20]. For the right-handed neutrino case ("type-I" seesaw), a comparison 
of the various calculations shows that there is no agreement on what is actually the result. 
This issue is also relevant for references using calculations in previous articles [19, 21]. Some 
of the discrepancies have a very significant impact on the predictions. Given both these 
experimental and theoretical situations, in this paper we recalculate the \i — )■ e conversion 
rate for the type-I seesaw model. We compare our results with previous calculations and 
determine the corresponding phenomenology it leads to. 

In the phenomenological analysis, we begin by comparing the conversion rate with the 
(j, — > e7 and [x —> eee rates; the latter process is very close to fj, — > e conversion in the sense that 
both transitions involve the same (local and non-local) form factors and consequently they will 
be shown to share common features in the asymptotic regimes; the former instead depends 
only on the dipole form factor for real photons, which also gives a non-local contribution to 
Li — > e conversion and ii — > eee decay. 

For a quasi-degenerate spectrum of right-handed neutrinos, a case which can naturally 
allow large rates, it turns out that the ratio of two Li-e processes (//—>• e in Titanium, li — > e in 
Aluminium, li — > cy, {i —> eee, ...) depends only on the right-handed neutrino mass scale [22], 
offering a quite interesting testing ground. As a function of that scale, we determine the 
Rn^e/ Br{jjL — > ej), R^ e /Br([i —> eee) ratios for Aluminium, Titanium, Lead and Gold 
nuclei, and also the R^t+J R^L^ e ratio. The observation of any two or more [i-e transitions, in 
agreement with the predicted ratios, would provide a strong evidence for such type-I seesaw 
scenarios. 

In this paper we also determine the range of right-handed neutrino masses and mixings the 
future [i — > e conversion experiments could reach. Although models with light seesaw scales 
may be unattractive theoretically as they tend to require large fine-tunings, it is worth to 
explore the full range of scales open-mindedly, as the results will apply to any singlet fermion 
in nature. It will be shown that the reach of planned [i — > e conversion experiments taken by 
themselves extends from 1000 TeV down to the MeV scale. Furthermore, for sterile species 
lighter than the electroweak scale, experimental bounds on the non-unitarity of the leptonic 
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mixing matrix in the electron and muon sectors [23, 24] become relevant and will be taken 
into account. For right-handed neutrinos lighter than ~ 2 GeV, the constraints on sterile- 
electron/muon coupling from K and D decays, as well as other constraints, superseed the 
sensitivity of fi — )■ e conversion processes. It will be shown that the planned [i — > e conversion 
experiments may become the main actor in detecting or setting constraints on sterile neutrino 
scenarios for mass scales from 1000 TeV down to 2 GeV. 

Sect. 2 describes the Lagrangian and couplings. Sect. 3 presents our calculation and 
deals with the general behaviour of the rates for the different processes considered, for any 
value of the right-handed neutrino masses; it also considers the particular cases of quasi- 
degenerate and hierarchical spectra. In Sect. 4 the regime of right-handed neutrino masses 
larger than the electroweak scale is discussed, both analytically and numerically, and the 
corresponding phenomenological analysis is performed; it also includes a comparison with 
previous computations in the literature. Sect. 5 focuses on the analogous analysis for heavy 
neutrino masses lighter than the electroweak scale, and in Sect. 6 we conclude. Exact analytical 
formulae and other items are included in the Appendix. 

2 Type-I seesaw model Lagrangian in mass eigenstate basis 



In the flavour basis, the Lagrangian of the type-I seesaw model with extra right-handed neu- 
trinos Nr reads 



C = C SM + iN R $N R 



N R Y N ftl L + ^N R MN R C + h.c. 



(2.1) 



where the flavour indices have been left implicit. Here Y^r denotes the neutrino Yukawa 
coupling to the Standard Model (SM) scalar boson (so-called "Higgs boson" for short), II 
denotes the left-handed lepton doublet, <f> is the higgs doublet and (j> = iT2<j)* with Ti the 
second Pauli matrix, 




(2.2) 



and v = 246 GeV. The charged lepton Yukawa coupling is assumed diagonal without loss of 
generality. We define the unitary (3 + k) x (3 + k) mixing matrix U through 




UP L n = UP L 



(2.3) 



V?, 

Ni 

where the vector n encodes all neutrino eigenstates in the mass basis, with Ni, i = 1 to k 
denoting the k extra physical heavy neutrinos, and Pl is the left-handed projector Pi = 
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(1 — 7s)/2- These mass eigenstates are Majorana fermions, n = n c , that is V{ = vf and 
Ni = Nf. In the mass eigenstate basis, the various gauge boson and scalar interactions read 
(for an arbitrary gauge, and with a, j3 denoting the three flavour indices, and i,j denoting the 
3 + k mass eigenstates) 

C W± = ^=W~ la^UaiPLTH + h.C. , (2.4) 

C Z = na^C^PLn, = ^-Z„ na 11 [C l3 P L - CtjP R ] rij , (2.5) 

= ~ Ji^ M $~ laUai (m ia PL ~ m ni P R )ni + h.c. , (2.6) 
^ = " *i C H ( m ^ p L ~ ™ nj Pr) fij , (2.7) 



igw ,3 _ 



_9w_ 
2M W 

gw 

4:M W 



[Cij (m ni P L - m nj P R ) - C*j (m ni P R - m nj P L )] rij , (2.8) 



h riiCij (m ni P L + m nj P R ) rij , (2.9) 
h rii [dj (m ni P L + m nj P R ) + C*. (m ni P R + m nj P L )] rij , (2.10) 



where g\y is the weak isospin coupling constant, c\y is the cosine of the weak mixing angle, 
and C and m n are (3 + k) x (3 + k) matrices defined as: 

3 

dj = ^2 U L U aj > m n = Diag(m„J = Biag(m l/1 ,m l/2 ,m U3 ,m Nl , ...,m Nk ) . (2.11) 

01=1 

The first three entries in m n denote the light neutrino masses, while the last k ones correspond 
to the heavy ones. Although the general results to be shown below will not assume any 
expansion in powers of Y^v/tun, for all plots we will neglect higher order contributions in 
this expansion. For this case it is worth to decompose the U matrix in blocks as 

which allows to recover the usual expressions valid up to 0[(Ynv, mi) /win] 2 corrections: 

Uuu = {1- € -)Upmns , U uN = Y N (M)~ 1 ^= , U Nv = -M- 1 Y N -j=U vv , U NN = , 

(2-13) 

where e = tj-Y^M _2 Y/V: e' = YjyY^M^ 1 , and Upmns denotes the customary lowest- 

order 3x3 leptonic mixing matrix appearing in charged currents. With these conventions, 
the neutrino mass matrix is given by: 

M V = -^Y%±Y N . (2.14) 
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3 fi to e conversion rates 



3.1 Calculation of the rates 

In the type-I seesaw framework, violation of charged lepton number arises at the one loop level. 
fj, to e conversion is induced by a series of gauge boson mediated diagrams given in Fig. 1. 
The various contributions to the process can be divided in those in which the momentum is 
transferred by the photon, by the Z boson or via two W bosons. The first two proceed via 
penguin diagrams, whereas the latter processes corresponds to a box diagram. Alike to the 
quark case, the internal fermions in the loop must have non-degenerate masses and non trivial 
mixings, in order to avoid a GIM cancellation. 

For a rigorous calculation of the rate it is necessary to separate the local contributions from 
the "extended" ones. This stems from the fact that extended contributions, unlike local ones, 
are sensitive to atomic electric field effects. The W and Z mediated diagrams are obviously all 
local. The 7 mediated diagrams contribute to both classes of transitions, extended and local. 
The /i — > e7 matrix element can be written as 



iM 



iegw 



2(4vr)2M^ 



4{q)u e (p') \F>} e {q 2 ^ - ^)P L - ia ^q u G» e {m e P L + m^P R )] u„(p) , (3.1) 



where q denotes the photon momentum, q = p — p' . The second term in this equation - 
mediated by the photon-lepton "dipole" G^ e coupling- is the only one contributing for an 
on-shell photon and is non local, whereas the "monopole" term F^ e is "local" (i.e. it only 
accounts for off-shell photon exchange and it involves 2 powers of the photon momentum in 
the numerator which compensate the long range 1/q 2 propagator of the photon between the 
lepton and nuclei lines [25]). One can therefore divide the effective Lagrangian relevant for 




(a) Photon Penguin Diagram (b) Z Penguin Diagram 



(c) Z Penguin Diagram 



W 



W 




(d) Box Diagram (e) Box Diagram 

Figure 1. The five classes of diagrams contributing to fj, to e conversion in the type-I seesaw model. 
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fj, — y e conversion in two parts 1 : 



= Sm a "* en > mF "' + w^£/ re ~ w * Yq + • (3 - 2) 

where F Xp = d x A p — d p A x denotes the electromagnetic field strength, and Fq C contains the 
contribution of the monopole F^ e term as well as that from the weak gauge-boson exchange 
diagrams. The coefficients Fg £ and G^ e encode all the dependence on the internal fermion 
masses and mixing angles. In Eq. (3.2), and all through the rest of the paper, the electron 
mass has been neglected. 

The effects of the nuclear form factors and of averaging over the atomic electric field can 
be taken into account in the way described in Ref. [26]. The nuclear information is encoded 
by D, and form factors whose values are shown in Table 1, taken from Ref. [26]. 
The final expression for the ratio of the fj, — >■ e conversion over the capture rate T cap t reads 



R 



2GW„m* 



(47r) 2 r, 



capt 



( 2 F pe + Ff ) + 4V™ (P pe + 2F% e ) + s 2 w G^Dj{2e) 2 (3.3) 



with, as usual, aw = 9w/^' K )i a = e2 /(^ 7r )) and Gf the Fermi constant. The expression for 
the dipole term G^f is given in the Appendix, Eq. (A. 25). The form factors Fq £ are in turn 
given by 



i> e = Q q s 2 w Ft? + Ff ( f - Q q s w j + \Fg« , (3.4) 

where q = u,d, Q q is the quark electric charge (Q u = 2/3, Qd = —1/3), X 3 is the weak isospin 
(Z^ = 1/2 , = —1/2) and sw is the sinus of the weak mixing angle. FJy £ , F^ e and Ff^ q 
denote form factors corresponding to photon-penguin diagrams, Z-penguin diagrams and box 
diagrams in Fig. 1. 

In the following we will use Eq. (3.3) for all numerical results. For low atomic number, 
aZ <C 1, and for comparison with other calculations, the results can be nevertheless simplified. 
In this case all the interactions can be considered point-like 2 (D/8e ~ V^) and the proton 
and neutron densities are approximately equal (V^ /Z ~ /(A — Z)), which gives 



R^e = -g^ »~^fF 2 p \Z (2*T + FT) + {A-Z) {Fj? + 2Ff) | 2 . (3.5) 



In this equation, A stands for the mass number, Z e ft denotes the effective atomic number, 
F p is a nuclear form factor related to through V^/V~Z ~ Z 2 ^F p a 3 ^ 2 /4ir, while F^ = 
(?) ^ Q<} s wGy e ■ The values of Z e ff, F p and T cap t that will be used in the phenomenological 
analysis below are given in Table. 1, taken from Refs. [26] and [27]. 



1 An interaction of the lepton current with the axial quark current is also generated but negligible compared 



with the dominant coherent fj, to e conversion, O(A) larger. 

2 The virtuality of the photon is set by the transferred momentum q 2 = 
considered point-like if this scale is small compared to the Bohr radius oc (Za) 
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Nucleus 


yip) 


y{n) 


D 


Z eff 


\F p {-ml)\ 


Tcapt (lO 6 ^ 1 ) 


£A1 


0.0161 


0.0173 


0.0362 


11.5 


0.64 


0.7054 


iTi 


0.0396 


0.0468 


0.0864 


17.6 


0.54 


2.59 


79 7 Au 


0.0974 


0.146 


0.189 


33.5 


0.16 


13.07 


jfPb 


0.0834 


0.128 


0.161 


34.0 


0.15 


13.45 



Table 1. Nuclear form factors and capture rates. 



Our complete analytical results for the amplitude of the various diagrams, together with 
their sum, are given in Eqs. (A.24-A.33) of the Appendix. The only approximations used in 
our calculation are to neglect: i) the electron mass compared to the muon mass; ii) higher 
orders in the external momentum over the W mass, as usual; iii) the value of the three light 
neutrino masses compared to the extra k heavier ones 3 . The first 2 approximations are 
accurate at a level better than O(10~ 4 ), that is to say better that what can be expected from 
higher loop contributions, and better than the uncertainties on the nuclear form factors. The 
third one becomes excellent as soon as the right-handed neutrino masses are a few orders of 
magnitude above the light neutrino masses. 

Finally, when comparing the strength of /i — > e conversion processes with the branching 
ratio for \i — > ej we will use for the latter the well-known result [28]- [32] 



Br(fi — > cy) 



256vr 2 T M 1 t 1 



(3.6) 



where r„ ~ 2.996 10 19 GeV denotes the total decay rate of the muon. When comparing 



instead with the strength of [i 
will be taken from Ref. [33], 



eee decay, the expression for the branching ratio Br{^ — > eee) 



Br(fjL — > eee) 



a, 



24576vr 3 Myy r M 



x <^ 2 



_ F ,eee + F ,e_ 2s 2 w{F ,e_ , e) 



+ IdslRe 



+ 324|G^ e | 



ln ^_ll 

ml 4 



+ 44|Ff -FJT\ 



484i?e [(Ff - )Gf *] 



(3.7) 



see the Appendix for the values of the various form factors. The comparison of this equation 
with Eqs. (3.3) and (3.4) illustrates that Br(^ — > eee) and [i — > e conversion are sensitive to 
different combinations of the same form factors. 



3 The latter approximation is not made in the first equality of Eqs. (A.24)-(A.33), only in the second equality 
of these equations. One could therefore use those results for any value of the right-handed neutrino masses. 



-7- 



3.2 Phenomenological framework 

The n — e transition rates expected in the type-I seesaw framework are clearly highly model 
dependent. For right-handed neutrino masses above the eV scale -as those considered in this 
work- for which the seesaw approximation U u n ~ Y^v/M in Eq. (2.13) holds, they scale 
in particular as the inverse of the right-handed neutrino masses at the fourth power, and 
contain four Yukawa couplings in the numerator, bringing each one a flavour dependence. 
Nevertheless it turns out that the models which can naturally give measurable rates are 
models which involve two or more quasi-degenerate right-handed neutrinos and for these 
models one can make remarkably clear predictions. The quasi-degenerate case is particularly 
natural, as it takes place for instance in scenarios in which lepton number (L) is approximately 
conserved [33-44]: the degeneracy is protected by the symmetry. This assumption allows a 
natural decoupling between large Yukawa couplings (inducing L-preserving large rates) and 
small Yukawa couplings (guaranteeing small neutrino masses) and results in viable scenarios, 
even for low seesaw scales. 

The remarkable predictions resulting for the quasi-degenerate case hold for the ratios 
of two rates where a same flavour transition occurs, for example 2?„_ >e /,Br(/i — > &y) or 
R^e/ Br(fj> — > eee). The point is simply that if the right-handed neutrinos are quasi degen- 
erate, only one right-handed neutrino mass scale is relevant in the rates and the dependence 
on the elements of the mixing matrix (which contain the Yukawa dependence) factorizes from 
the mass dependence. For instance, for ~ mjv 2 ~ ... = m^, the rates in Eq. (3.5) for 

fj, —> e conversion in light nuclei and Eq. (3.6) for [i — > ej [the generalization to heavy nuclei 
is straightforward from Eqs. (3.3) and (3.5)] take the factorized form 

K ^ e S^T capt Z ** 

k 

(A + Z)F u (x N ) + (2A-Z)F d (x N )} 2 \Y^U eNi U* N ^ . (3.8) 

i 

Br^ey) ^ |E^^f > (3-9) 

where G 1 {x^), F u (xn) and F^xn) are functions of x^ = m 2 N /My / given in Appendix, 
Eqs. (A. 9), (A. 42) and (A. 43). In consequence, the Yukawa coupling (flavour) dependence 
drops in the ratio of both rates, leaving only a dependence on the mass tun of the heavy 
neutrinos [22]: 



(3.10) 

Note that the cancellation of the Yukawa coupling dependence in the ratio is valid at the 
dominant order in the mixing angle expansion, that is to say that in Eqs. (3.8)-(3.10) we 
neglected terms which involve four insertions of light-heavy mixing in the amplitudes, see 



(A + Z) F u (x N ) + (2A - Z) F d {x N ) 



1 2 



G 7 (x N ) 
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Appendix. This approximation is justified for the range of low right-handed neutrino masses 
we contemplate, in view of the experimental constraints on the mixing, whereas for the large 
mass regime it relies on the perturbativity of the Yukawas. Note also that, for a given nucleus, 
the ratio may vanish for a particular value of mjvj see Sec. 4 below. 

The predictions on ratios of rates, such as that in Eq. (3.10), hold as well for the case where 
only one right-handed neutrino would dominate the rates, because in this case the Yukawa 
dependence also drops, leaving the same dependence on a unique mass. This situation is quite 
generic of scenarios where the right-handed neutrino mass spectrum is hierarchical, since in the 
large > My/ (low tum < My/) mass regime all but the lightest (heaviest) Ni contributions 
can be in general neglected, see Eqs. (4.6), (4.7), (5.6) and (5.7). However, the price to obtain 
measurable flavour-changing rates with hierarchical spectra is to induce unacceptably large 
neutrino masses, disregarding eventual large fine-tuning between the various parameters 4 
(which anyway would be destabilized by unacceptably large radiative contributions to neutrino 
masses [46]). Moreover, the heavy right-handed neutrino contribution to neutrinoless double 
beta decay can be in contradiction with the current bound [47-49], especially for low right- 
handed neutrino masses below a few tens of GeV. In consequence, in the rest of the paper we 
will consider only the quasi-degenerate spectrum of right-handed neutrinos. 

Although the full range of right-handed neutrino masses has been analyzed in this work, in 
the next sections we discuss in detail the maximum and minimum right-handed neutrino mass 
scales that future sensitivities may reach, as well as the sensitivity reach for the charged-current 
mixing of steriles with the electrons and muon sector of the SM. We will denote by "large" 
mass regime that in which the right-handed neutrino scale is larger than the electroweak scale, 
mjv > My/ , while the "low" mass regime will be that in which the heavy right-handed neutrinos 
are lighter than My/ (although always much larger than the usual three light neutrinos). 

4 Large Mass Regime (uin > My/) 

Let us consider first in detail the regime of singlet fermion masses larger than the electroweak 
scale. Before reporting on the analytical results for G^ e and F^ d , it is interesting to discuss 
the expected behaviour of the transition rates mediated by extra neutrino species much heav- 
ier than My/. In the limit of infinite right-handed neutrino masses, the low-energy theory 
is renormalizable (e.g. the Standard Model with massless left-handed neutrinos) and, in con- 
sequence, the extra singlet degrees of freedom introduced must decouple, leaving no impact 
on low-energy observables, see for example [50]. We have checked explicitly that indeed all 
rates mediated by right-handed neutrinos fulfill this condition. As an illustration, consider 
the Z-mediated contribution to Fq £ in Eq. (3.3). Although the loop integral exhibits a loga- 

4 For instance in pseudo-Dirac models as in Refs. [44, 45], if one increases the mass splitting, one gets 
too large neutrino masses, unless one would have a cancellation between the different types of neutrino mass 
contributions, i.e. from the Yukawa couplings whose simultaneous presence breaks lepton number and from the 
two right-handed neutrino mass matrix entries which break lepton number (the one contributing to neutrino 
masses at tree level and and the one contributing radiatively, see for example Eqs. (3)-(9) of Ref. [45]). 
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rithmic m? N /Myy dependence, and furthermore longitudinal Ws in the loop provide positive 
powers of mjv, the 1/M (implying 1/mjv) dependence of the elements of the mixing matrix 
-see Eq. (2.13)- ensures a total rate scaling as 1 / rri^log 2 {rn 2 N / My^) for S> M\y. In other 
words, a rate that vanishes in the decoupling limit as it should. See also the next subsection. 

4.1 Analytical results 

Expanding Eqs. (A.24-A.33) to lowest order in inverse powers of = mr N ./Myy, the terms 
relevant for the ratio of the fj, — >■ e conversion rate to the capture rate, Eq. (3.3), read 

k 

K e = £ u eNi u; Ni F u ( XNi ) , f u {x) 
i=i 

k 

K e = £ u eNi u; Ni F d { XNi ) , F d { x ) : 
i=i 

k 

= £ UeN.U^G^XN,) , G 7 (x) 
1=1 

whereas for the particular case of light nuclei in Eq. (3.5), it results 

k 

K e = £ U eNi U; Ni F u (x Ni ) , F u (x) 

i=l 
k 

F T = E U eN>U; Ni F d ( XN J , F d ( X ) 
i=l 

Comparison with the literature 

Although the first two items below hold for the exact formulae and not only for the limit of 
large singlet fermion masses, we gather here together the result of comparing our formulae with 
those in previous literature (most often given in the approximation of low atomic number). A 
number of comments can be made: 

• Sign of the G 7 contribution. Our results agree with the sign indicated in Refs. [11, 15, 16], 
and is opposite to that in Refs. [12, 13]. Note that G 7 contributes both to fj, — > e 
conversion and to \i — > ej decay -see Eq. (3.1), and the loop integrals involved can be 
related with analogous amplitudes in the quark sector, once the internal quark charge 
is switched off: this allows to check that our results are consistent with those for K 
transitions [51] and for b —> sl + l~ decay [52]. 

• Box diagrams. Refs. [11-13] also exhibit differences in the relative size and/or sign of the 
crossed and not crossed box contributions, that is, in the amplitudes resulting from the 
last two diagrams in Fig. 1. As shown in the Appendix, we obtain a —4 factor between 
these contributions, in agreement with for instance the results for AS" = 2 or AB = 2 
processes in Ref. [53]. Also, we obtain half the contribution considered in Ref. [21] for 
all box diagrams. 



2 a 16 log (a) - 37 3 + 31og(x) 
3 Sw 12 8 (4 ' 1] 



1 2 16 log (a) - 37 3- 31og(x) , 
'3 Sw 12 8^ ' ' (4 ' 2) 



(4.3) 



2 2 16 log (a?) -31 3 + 31og(x) 
3 Sw 12 8^ ' ' (4 ' 4) 



1 2 16 log (x) -31 3-31og(aQ 
3 Sw 12 8^ 



.(4.5) 



• Constant terms in the Fz et form factors. For right-handed masses low enough to 
result in observable n — )■ e rates, those terms are numerically competitive with logarith- 
mic ones and cannot be neglected, see Eqs. (4.4)-(4.5). We get different results for them, 
though, than in Ref. [16], where some of them have been neglected. On the contrary, 
our results agree with the full expressions for those form factors given in in Ref. [33] for 
the computation of the /i — > eee rate, to which they also contribute. Furthermore, the 
fi — > e conversion rate has been also calculated in the framework of the supersymmetric 
type-I seesaw model in Ref. [15]; we checked that its non-supersymmetric limit results 
in logarithmic terms which agree with ours, but there too the constant terms have been 
neglected 5 . 

• Decoupling limit. As already mentioned above, in the limit in which the right-handed 
neutrino masses go to infinity, Eqs. (4.4)-(4.5) show that these states decouple in all 
rates considered, as they should [50]. This is in disagreement with the non-decoupling 
behaviour obtained in Ref. [19] for the fi — > e conversion rate in the type-I seesaw 
model 6 . Note that the heavy- light mixing angles are suppressed by one power of the 
right-handed neutrino masses, i.e. U ~ Yvjm^. In summary, in the large mass regime 
considered here, xn = m%/Myy 3> 1, the leading term in the three type of rates under 
discussion scales as 

r ~ (logXAr) 2 /x^r , for /i — > eee and //—deconversion, (4.6) 

r~l/a&, for/i^e7. (4.7) 

Finally, note that with respect to the results in Ref. [14] we get different coefficients for 
several logarithmic and constant terms. 

4.2 Ratios of rates involving one same flavour transition 

The formulae in Eq. (3.3), (3.6) and (3.7) allow to compare the relative strength of fi — > e 
conversion to the fi — > e*y branching ratio, R^Zf^, and to the /i — > eee branching ratio, R^Z^ eee . 
As explained above these ratios depend only on the right-handed neutrino mass scale thn (at 
leading order in Yv/M). The results are illustrated in Fig. 2 as a function of win for various 
nuclei, in the large m^r regime 7 . One distinctive feature of these ratios is that they vanish for 
some value of tun, when the [i — > e conversion rate, Eq. (3.3), vanishes. This peculiar feature 
is due to the up quark Fu £ and down quark F^ e contributions in Eq. (3.3) having opposite 

5 Numerically, to neglect the constant terms for these form factors leads for example to a Titanium con- 
version rate vanishing for mjv — 2.9 TeV instead of 4.7 TeV, which translates into rates which can differ by 
several orders of magnitude in the few TeV range. 

6 The analytic expressions used in Ref. [19] are the same ones that are valid for the case of a fourth 
generation of quarks and leptons [54] (that is, with an active extra neutrino instead of a singlet right-handed 
neutrino, which cannot decouple as the remaining low-energy theory would not be renormalizable) . 

7 For Mjv w 100 GeV, we get R^ e w 6.7 • Br(n -» erf) and Br{^i -> eee) w 0.033 • Br{^i -» ey), which 
differs with Eq. (3.21) of Ref. [16] by two orders of magnitude. 
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m N [GeV] m N [GeV] 



Figure 2. R^Xe-y = R^->e/Br(n — > ej) (left panel) and R^Xtee = R^e/Br(fi —> eee) (right panel) 
as a function of the right-handed neutrino mass scale mjv, for /i — > e conversion in various nuclei. 



signs, as an outcome of their different charge and weak isospin. The precise value where the 
fj, — > e conversion rate vanishes is nuclei-dependent and given by 



which shows that small variations on the nuclear form factors may result in sizeable variations 

on the value of m? N , which is thus sensitive to the nuclear physics uncertainties. The uncer- 
o 

tainty in the ratio /V^ translates, for instance for 5-10% variations, into O(TeV) shifts 
on the value of the right-handed neutrino mass at which the conversion rate vanishes. With 
the form factor values given in Table 1, the rate vanishes for mass values typically in the 2-7 
TeV range, respectively 6.4, 4.7, 2.5 and 2.4 TeV for Al, Ti, Au and Pb, as Fig. 2 shows 8 . 

For degenerate right-handed neutrinos and light nuclei, aZ <C 1, Eq. (3.8) is a good 
approximation which allows to rephrase the vanishing condition as 

F u (2A - Z) 

T d ~~ (A + z) ■ (4 ' 9) 

Due to the logarithmic behaviour of F u /Fd, a small variation of {2A — Z) / (A + Z) results 
in a sizeable variation of the value of tun for which R^ e vanishes. The atomic ratio on 
the right-hand side of Eq. (4.9) takes the value —1.05 for Ti , —1.02 for Al, and —1.15 for 
Au and Pb. For illustrative purposes, Fig. 3 shows the value of F u /Fd as a function of m^, 
together with the value of (2^4 — Z)/(A + Z) for each nucleus. Furthermore, performing in 



Note that a plot of the same ratio is displayed in Ref. [19], with quite different results, in particular 
vanishing rates for much lower mjv values, see footnote 6 above. 




°' 9 2000 4000 6000 8000 10000 
m N [GeV] 

Figure 3. \F U /Fd\ (that is, —F u /Fd) as a function of the right-handed neutrino mass scale (solid 
line). The dashed lines give the values of {2A — Z) / (A + Z) for the various nuclei. Crossing points 
give the values of where the fi — > e conversion rate J?u_> e vanishes. 



this approximation an expansion in inverse powers of xn 
and (4.5), allows to rewrite the mass value at which 

%A-Z) + 



m 



/Myy, as given in Eqs. (4.4) 



N 

vanishes as 



m N 



M w exp 



9 , 314, 
8 ' 12 



\(A-z) + 



4 S 2 



z 



z 



(4.10) 



Fig. 4 depicts two other ratios for /j, — > e transitions: Ril^. e / Rwi^. e and Br(fi — > e-y)/Br(/j — ; 
eee) (the latter one from Ref. [22]). Sweeping over increasing mjy values, the first ratio first 
does, and later goes to infinity when R^L> e vanishes. The second ra- 



vanishes when RzL^ e 
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Figure 4. R^% e /R^ e (black) and R^Jele — Br(fJ, — > e^)jBr{[i — » eee) (red) as a function of ■ 
The horizontal dashed line show the large asymptotic value of RTL ye /R^L > . e . 
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tio, Br(pt — > cy)/Br(n — > eee), is monotonous in mjv, hence it does not display any tun 
degeneracy. 

The ensemble of the results above imply that, from the experimental determination of 
two fj, — e transition processes, and up to discrete degeneracies, it is possible to determine the 
scale vci jy of the generic framework considered. That pair of processes could be any two among 
the four processes which will be probed with improved sensitivity in near future: /x —> ej, 
fx -> eee, Rf_, e and R^ e . 

To lift possible degeneracies a third measurement may need to be considered. As an 
example, assume that from the MEG and COMET experiments, R£ _> e .y is measured to be 
~ 0.1: the values mjy ~ 2.5 TeV or mjy ~ 16.5 TeV would then be singled out, see Fig. 2. To 
lift this degeneracy the observation of a third /i — > e transition process would be necessary: 
for instance Rj^l^ e at PRISM or fi — > eee at fj,3e — PSI [55]. Alternatively, the measurement of 
two rates might be incompatible with the upper bound or measurement of a third one, which 
would rule out the scenario 9 . Similarly the measurement of a single rate, together with the 
upper bound or measurement of another one, could exclude this scenario for ranges of m^j 
values (eventually excluding the whole mass range). 
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Present bounds and future sensitivity to (b) Present bounds and future sensitivity to 
Yukawas as a function of the heavy-light mixing as a function of the 
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Figure 5. Bounds on | Y^l ^Vi,»^iV* 1 an d U e NiU* N .\ for scenarios characterized by one right- 
handed neutrino mass scale. The solid lines are obtained from present experimental upper bounds: 
from Eq. (1.3) and Br{{± ej) < 2.4 • KT 12 [56], Br{{i ->■ eee) < KT 12 [57]. The dashed lines 
are obtained from the expected experimental sensitivities: from Eqs. (1.1)-(1.2) and Br(/i — > cy) < 
10- 14 [58], Br{fi -> eee) < lfr 16 [55]. 



9 Note also that, analogously, the measurement of r — > lj decay and of r — > II' I' decay would also allow to 
determine the mjv scale [22]. That determination could be compared with the fj, e results above, to rule 
out or further confirm this scenario. 
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4.3 Maximum seesaw scales that future experiments could probe 

Fig. 5 shows the lower bounds resulting for the Yukawa couplings and mixing parameters, if 
the various rates are required to be large enough to be observed in planned experiments. It 
also shows the upper bounds which hold today on these quantities from the non-observation 
of these processes. This figure illustrates well the impact of future [i — > e conversion measure- 
ments/bounds, as they will become increasingly dominant in exploring flavour physics in the 
\i — e charged lepton sector. Values of the Yukawa couplings as low as 10 _1 , 10 -3 and 10 -4 
could be probed, for tun = 100 TeV, thn = 1 TeV and = 100 GeV, respectively, with 
Titanium experiments being the most sensitive. If the Yukawa couplings are required to lie in 
the perturbative regime, i.e. that each Yukawa coupling is smaller than ~ V^tt, the bounds 
of Fig. 5 can be rephrased as upper bounds on the tun scale: 

m w < 6000 TeV- (-^—) 4 , (4.11) 
m N < 1000 TeV- — =- 4 , (4.12) 

io- 14 \ \ 



^< 300IteV-( j , (4.13) 



IO" 16 \\ 



m N < 1000TeV- -) . (4.14) 

\Br([i — > eee)/ 

Imposing instead that the Yukawa couplings should be smaller than unity would lead to 
bounds smaller by about a factor of 3. Overall, this exercise shows that future experiments 
may in principle probe the type-I seesaw model beyond the ~ 1000 TeV scale. 

5 Low Mass Regime (m^v < Mw) 

This section focuses on the low mass region (tbat < Mw), to asset the discovery potential to 
singlet fermions expected from future /i — > e conversion, [i — > e^f and [x — > 3e experiments. 

5.1 Analytical results 

In the low mass regime, expanding in powers of the small parameter xjy = m? N /M w <C 1, 
the leading terms of the different form factors relevant for /i — > e conversion in an arbitrary 
nucleus -see Eq. (3.3)- are given by 



i?^ ST TT rr * £, / x „/x (1 2 41og(a:)+6 3 + 61og(x) ^ 

K = Z^ UeN i U nNi F u( XN i-> ' F ^ x > = [s s W 4 + g J x i I 5 - 1 ) 

k 

*T = E U eNl U; Ni F d (x Nt ) , F d {x) = (-^ 41 ° g( 4 )+6 +|)^, (5-2) 



i=l 
k 



^ = E UeN^G^XN,) , G 7 (x) = - , (5.3) 



i=l 
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whereas for conversion in light nuclei, Eq. (3.5), they take the form 



FT = t v^Aim) • *(*> = (|4,U!lM±I + S±«afia) x , (5.4) 

i=l \ / 

*T = £ tU.^F^JVi) , W = (-^ 4l0g( 4 X) + 7 + I) x • (5.5) 

i=l ^ ' 

As a consequence, for xn = m %/M^y <C 1 the leading terms in the transition rates vanish 
as 10 

T ~ x%(logXN) 2 , for fj, — > eee and /i — > e conversion, (5.6) 
T ~ x 2 N , for [i — > e7 . (5.7) 

This is in contrast with the leading behaviour found for the large mass regime xn S> 1, 
Eq. (4.6), with the scaling law being inversely proportional to xf^, ensuring decoupling. Note 
that to get Eqs. (5.6)-(5.7) we assumed fixed U v n mixing parameters, as it is customary to 
express constraints on sterile neutrino models in terms of these mixing parameters 11 . 

5.2 Minimum seesaw scales that future experiments could probe 

Taking /j, — )■ e conversion experiments by themselves, the sensitivity to low singlet fermion 
masses is outstanding. This is illustrated in Fig. 6a (dashed lines), which shows a sensitivity 
down to 2 MeV for Titanium for large mixing. Nevertheless, a series of existing experimental 
bounds curtails the expected impact, to wit: 

• Unitarity bounds. The mixing matrix elements Ugw entering the rates are constrained by 
the bounds on the unitarity of the leptonic mixing matrix [23, 24]. The relevant bounds 
here are those holding for mjy < M\y, in which the strong constraint from fi — > cy 
present for masses above Myy is lost due to the restoration of the GIM mechanism under 
the M\y scale. In the region ~ 100 MeV< < M\y the bounds are dominated by the 
constraints on f3 decay (lost under the GeV regime), kaon decays, universality constraints 
from r and tt decays, and tree-level [i decays. For very light N{, under the pion and 
muon masses, the constraints from the decays of the latter are lost, since all eigenstates 
are then available in the decay, so that unitarity is recovered. Nevertheless, below 100 
MeV the absence of zero distance effects in short baseline oscillation experiments (such 

10 In the low mass regime, the amplitude for /i — > ej is analogous to that for b — > sj, while those for fj, — > e 
conversion and [i — > 3e exhibit only a GIM cancelation quadratic in the light masses instead of the logarithmic 
one for quark transitions such as 6 — > se + e~, which are proportional to fermion electric charges inside the 
loop. 

11 For small sterile masses above the eV, the seesaw approximation U v n oc Y^v/itin (<S 1) in Eq. (2.13) 
still holds, and an extra factor has to be added in Eqs. (5.6)-(5.7), so that the rates have a logarithmic 
or constant dependence on mjv. Below this range one enters the Dirac-dominated regime, where the mixing 
angles become free parameters, and the asymptotic behaviour is that in Eqs. (5.6)-(5.7). 
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m N (GtV) 

(a) Sensitivity reach of present and future 
experiments in the mixing-mass 
parameter space 

Figure 6. The left panel shows the | J2i U e NiU* N . | versus mass sensitivity regions for present (contin- 
uous curves) and future (dashed curves) e — /i flavour experiments. Black, red, green and blue curves 
result from Br(/i — > ej), Br{[i — > eee), R^L^e and R^L,. e , respectively. The regions already excluded 
by non-unitarity limits, ir and K peak searches, 7r, K, D, Z decay searches, BBN, SN1987A and LHC 
collider searches (dotted lines) are also indicated. Shaded areas signal the regions already excluded 
experimentally. The right panel shows the maximum allowed flavour changing rates compatible with 
the bounds of the first panel. The horizontal lines give the present (solid [7, 56, 59]) and future 
(dashed [2-4, 55, 58]) sensitivities of the different experiments. 

as KARMEN [60] and NOMAD [61] ) also sets constraints on the mixing elements [23]. 
Stronger constraints follow nevertheless in that region, mainly from "peak" and "decay" 
search experiments. 

• Peak experiments explore the direct production of light (< mx, Tn^) extra singlet fermions 
in two-body (INi) particle decays of light mesons. From pion [62] and kaon decays [63, 
64], the absence of a monochromatic line -or peak- in the charged lepton energy spec- 
trum at (w-Xjt + wJ -m^)/2mx,7r 12 , excludes at present the 30 MeV< mjy < 400 MeV 
region. Decay searches provide even stronger constraints. 

• Decay experiments including more than 2 particles in the final state look for the effects of 
the production and decay of massive neutrinos. The relevant processes for constraining 
J2i U e N,U* Ni are K,ir,D -)■ INi -> tl'vyt with £, I" = e, fi. Their non observation 
sets very strong constraints in the range 1 MeV< m^r < 2 GeV [65-68]. Similarly 
searches for a Z — > NiU decay sets interesting constraints below the Mz mass [69]. For 



m N [GeV] 

(b) Maximum rates allowed by the various 
existing constraints 



12 In the rest frame of the decaying meson. 



a more detailed discussion of both decay and peak experimental bounds see Ref. [70], 
and also Refs. [71] and [72]. 



• Supernovae and BBN limits. Upper and lower bounds on the mixing can be obtained 
from supernova SN1987A data [73], since the known duration of the blast would be 
modified if the right-handed neutrinos are produced in the core and escape carrying 
energy away. BBN limits have also been explored for the very low mass region [72, 74- 
76]. 

• LHC search for the decay of a 125 GeV Higgs boson via the process u N — )• v uli provides 
new bounds on the seesaw parameters [77-80], which can be translated into bounds on 
rare processes. The branching ratio of the Higgs boson decaying into Majorana neutrinos 
enters the implementation of these bounds, which in the general seesaw here considered 
reads 



Br{h -> vN) = °" {\U eNi \ 2 + |£W 2 + \U m \>) m h m% 1 - ^ ■ , 

° IV1 W L hi \ m h ) 

(5.8) 

where the decay into all three light neutrinos Vj and all heavy ones Ni is considered. r^°* 
denotes the total decay width of the Higgs boson including the SM channels plus those 
producing the extra heavy neutrinos. The different decay widths of the heavy neutrinos 
also enter the analysis, mainly all N — > t>l + d,~ channels. 

Following Ref. [70] (see also Refs. [71, 72, 76, 81]), Fig. 6a provides an approximate drawing 
of the regions excluded by these constraints (depicted as shaded areas), in particular non- 
unitarity bounds, kaon and pion peak searches, kaon decay searches by the PS191 experiment, 
D meson decay searches by the CHARM [67] and NuTeV [68] experiments 13 , Z decay searches 
by the Delphi experiment [69], and Higgs decays data from LHC. The exclusion lines for the 
mixing parameter \T,iU e N i U* N . | are valid at 90% CL, for the present (continuous curves) and 
future (dashed curves) reach of the three types of measurements under discussion: /i — > ej, 
[i —> eee and jj, — > e conversion. Note the impressive sensitivity expected from [i — > e 
conversion in Titanium taken by itself, reaching seesaw masses down to 2 MeV. Nevertheless, 
the very stringent PS191, CHARM and NuTeV bounds on decay searches into N{ determine 
a lower bound on the mass at which conversion experiments would be competitive: 2 GeV. 

Fig. 6b shows the maximum rates allowed by the bounds on mixing depicted in Fig. 6a 
(shaded areas), as a function of the mass m^, at the 90% CL. Horizontal lines in Fig. 6b 
indicate the present (continuous) and future (dashed) sensitivities of the different experiments: 



2 x 2 



13 The bounds from meson decay searches for masses below 2 GeV depicted in Fig. 6 are those provided 
by the experimental collaborations which assume decaying Dirac neutrinos; for Majorana neutrinos they may 
differ by factors ~ \/2 or less, depending on whether the decay channel is or is not self-conjugate [65]. The 
inclusion of those factors would require a reanalysis of the experimental data for masses lighter than 2 GeV, 
which we refrain from attempting here as that region turns out to be out of reach for the /j, — > e conversion 
experiments under discussion. 



their intersection with the corresponding rate (depicted with the same color) determines the 
lowest mass for which the new experiments would improve present bounds. 

About the Higgs boson decay constraints, Fig. 6a depicts a model-independent bound, 
except for the assumption of degenerate heavy neutrinos (or hierarchical ones, e.g. dominated 
by only one heavy species). This bound is obtained [79] by translating the absence of an excess 
over the SM expectation for the channel h — > vvt r t~ into an upper bound on Br(/i — > vN) < 
0.4 [82] in Eq. 5.8, and employing then the general inequality lEjC/eA^t^y. | < ^i. a \U a Ni\ 2 
. Figs. 6a and 7 illustrate that it turns out to be less stringent than for example the 
present MEG bound on Br(fi —> e^f). In the future, if LHC can reach a ~ 1% sensitivity on 
Br(/i — > vN), the constraint would be comparable to the present MEG one. 

More stringent bounds follow in concrete realizations of the seesaw. For instance, the 
scenario with just two right-handed neutrinos added to the SM and approximate L conser- 
vation in Ref. [44, 45] (see also Refs. [38-40, 43]) is very predictive. The mixing-dependence 
of Br(/i —7- uN) -the term between parenthesis in Eq. 5.8- reduces in this case to an overall 
scale dependence, which can then be bounded from Fig. 6a. Furthermore, in these scenar- 
ios the mixing elements l^aivj are not arbitrary, but explicit functions of a few observable 
quantities such as the measured light neutrino mass differences and mixing angles, the over- 
all scale and the Dirac {5) and Majorana (a) phases. This allows to express bounds on 
\U e NiU* N . | from Higgs decay as a function of the values of those CP phases. This inter- 
esting fact is illustrated in Fig. 7, whose bands depict the maximum and minimum bounds 
obtained by varying those phases. The values of the neutrino parameters used in this figure 
are Am^ = (7.59 ± 0.2) • 10~ 5 eV 2 , (Am 2 ^ = (2.36 ± 0.2) • 10~ 3 eV 2 for the case of neutrino 
inverted hierarchy (IH), lAm 2 ^! = (2.46 ± 0.12) • 10 -3 eV 2 for the case of neutrino normal 
hierarchy (NH), 6 l2 = 34.4 ± 1 degrees, 6» 23 = 42.8 j^g degrees, and sin 2 26 13 = 0.092 ± 0.017 
[83, 84]. For NH, the maximum and minimum boundary lines of the bands correspond approx- 
imately to (5, a = 0,7r/2) and (5, a = 0, — vr/2), respectively, while the equivalent values for 
the case of IH are (<5, a = 37r/2,7r) and (5, a = 0, — vr/4). The values allowed never reach the 
absolute bound, also depicted. On the other side, note that the plot for IH allows some points 
in which the fi — e mixing tends to vanish: this is expected and well-known (first pointed out 
in Refs. [85, 86]), because in this class of models some entries of the Yukawa couplings may 
vanish (and therefore the entry in the mixing U a ^ i would vanish too) for certain values of the 
CP phases and light neutrino mixing and mass differences within their 90% C.L. region [22]. 
As an example of the impact of the analysis for this type of models, with the present limit 
on Br(/i — > vN), if MEG observes a signal of order the present sensitivity, only particular 
ranges of the phases would be allowed in this model. In any case, for the /U-e channel it is not 
expected that LHC data will allow to approach, even from far, the sensitivities one expects in 
the new generation of /i — > e conversion experiments, Eqs. (1.1) and (1.2). 

14 In Fig. 6a this bound is plotted down to 20 GeV, a point below which the typical cuts in the invariant 
mass of the lepton pair would remove the events coming from higgs decay via the right-handed neutrino [79] . 
The equivalent bounds shown in [79, 80], which we used, do not extend further down than 50 GeV: in the 
region 20 — 50 GeV what Fig. 6a depicts is an extrapolation. 




Figure 7. Bounds on lUiU^U^. | from LHC Higgs decay data and comparison with flavour violation 
searches. The contraints from LHC data are illustrated for a sensitivity of Br(h —> vN) < 0.4, and 
a future one improved to the % level [80]. For these values, isolated green curves are the absolute 
bounds for a generic seesaw model. Bands in orange and red show instead the variation with the 
unknown values of the Dirac and Majorana CP phases, for the approximately L conserving scenario 
in Ref. [44], for normal (lelft panel) and inverted (right panel) hierarchy. The present (black) and 
future (dashed black) MEG sensitivities and the expected one for conversion in Titanium (blue), in 
Aluminium (green) and /i — > eee (red) are shown for comparison. 



Finally, Fig. 8 depicts the ratio of fj, — > e conversion rate to Br{[i — > e^f) and to Br(fj, —> 
eee), in the low mass region of the type-I seesaw scenario considered. These ratios are now 
always larger than one, and they do not vanish for any value of the singlet fermion mass, in 
contrast with the behaviour in the large mass regime discussed in Sect. 4. With diminishing 
mass, the ratio B^X% ee becomes constant, while the ratio R^^%y grows logarithmically, see 
Eqs. (5.1)-(5.5). The dotted section of the curves indicates the mass region which is out of 
observability reach for the planned experiments under discussion. 

Note that in Fig. 8 the ratios to \x — > e"f and fi — > eee for Titanium and Lead fall on 
top of each other. Numerically the agreement holds at the % level. This is in fact due to the 
conversion ratio R^ e in both nuclei being approximately equal, at the percent level (modulo 
possible nuclear physics uncertainties). The ratios entering this quantity vary from atom to 
atom solely through the variables / \jT cap t , V^ 1 ' / \jT cap t , D / \fT cap t, see Eq. (3.3), which 
differ by about 10% for Lead versus Titanium, and similarly for Gold for instance, see Table 1. 
But the combination in which these variables enter in Lead and Titanium, via the form factors 

,Gf^ e in Eq. 3.3, happens to be very close numerically. The coincidence at the % level of 
the fj, — > e conversion ratios R^ e in Titanium and Lead in all the mass range under 50 GeV is 
a prediction of the quasi-degenerate right-handed neutrino scenario for very low seesaw scales. 

It is worth to remark that the sensitivities and exclusion plots obtained apply not only 
to the (quasi-degenerate) type-I seesaw scenario object of this work, but to any BSM renor- 
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(a) Ratio of the conversion rate to fi — > (b) Ratio of the conversion rate to fj, — > eee 



Figure 8. fi — > e conversion rate over Br(fi — > &y) (left panel) and over Br(fi — > eee) (right panel) as 
a function of the right-handed neutrino mass scale mjv, in the light regime (m^v < M\y) for various 
nuclei. Curves become dotted lines at the limit of observability by future experiments. 

malizable theory containing Dirac or pseudo-Dirac singlet fermions, which mix with SM light 
leptons with strength U^. 



6 Conclusions 

Future experiments aiming to detect fi —> e conversion in atomic nuclei are especially promising 
for the discovery of flavour violation in charged-lepton transitions. For various nuclei, we have 
calculated the — > e rate which holds in the framework of the type-I seesaw scenario of 
neutrino masses. We discussed the phenomenological impact and the reach of present and 
planned experiments in terms of the seesaw scale and mixing. 

Analytically, the fj, — > e conversion rate includes form factors with and without a log- 
arithmic dependence on the heavy singlet fermion masses. For the former (unlike for the 
latter) our results agree with those in Ref. [16] (and in Ref. [15] provided one takes the non- 
supersymmetric limit of the results quoted there). We basically disagree at different levels 
with all other calculations we have found in the literature, be it for the constant or the loga- 
rithmic terms. The constant terms in the form factors turn out to be numerically competitive 
with the logarithmic ones and must be taken into account. 

The computation performed is valid for all values of the heavy right-handed neutrino 
masses above the MeV scale (that is, sizeably larger than the masses of the three known 
neutrinos), and the complete expressions have been used to obtain the numerical results. 
Furthermore, for illustrative purposes, the leading behaviour in the regime of large (> Myy) 
and of small (< Myy) right-handed neutrino mass scales has been explicitly discussed. 
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The various possible ratios of rates involving the same charged /x — e flavour transition, 
have been determined. When a single heavy scale dominates, they depend exclusively on that 
unique scale. The results are illustrated in Figs. 2 and 4 for the large mass regime, and in 
Fig. 7 for the low mass regime. Fig. 9 provides a summary of these results by merging these 
figures together. It is useful for comparison purposes and illustration of degeneracies within 
one same experiment and across experiments, especially if a positive signal is observed. These 
results apply both in the limit of non-degenerate heavy neutrino masses and in the hierarchical 
limit. We focused the discussion on the quasi-degenerate scenario, as it is the most natural 
one to allow for observable rates. 
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Figure 9. Ratio of the /i — > e conversion rate in Al (green lines), Ti (blue lines) and Au (red lines) 
to Br(fi — > e"f) (solid lines) and to i?r(/x — > eee) (dashed lines) for the entire range in mass tun here 
considered. Lines are dotted when they require, for /i — > ej and /i — > eee, a sensitivity better than the 
one expected at planned experiments. 

Since the various ratios exhibit a different mass dependence, they offer a particularly neat 
way to confirm/exclude possible scenarios. In the case of an agreement between theory and 
measurements, they would allow a determination of the right-handed neutrino mass scale. 
For right-handed neutrinos in the large mass regime, depending on the ratio which would 
be measured, one may have a degeneracy in the right-handed neutrino mass scale it could 
correspond to. This can be resolved by the measurement of another ratio. Furthermore, we 
find that the values of the heavy neutrino masses for which the ratios of \i — > e conversion to 
ji decays may vanish are around 2 — 7 TeV. 



We have explored the sensitivity to the singlet fermion mass scale as a function of the 
mixing with the muon and electron sector of the SM. The maximum scale that future [i — > e 
conversion experiments could probe is above the 1000 TeV scale. This sensitivity extends to 
very low masses, down to ~ 2 MeV for Titanium. When the stringent bounds coming from 
massive sterile neutrino searches (from the unitarity limits of the Upmns mixing matrix, from 
7r, K and D meson decay searches, and other constraints) are taken into account, the expected 
impact in constraining the low mass regime is somewhat reduced. We have discussed as well 
the LHC sensitivity from Higgs decay via heavy neutrinos, both in a model-independent way 
and for the case of seesaw models with approximate lepton number symmetry, comparing it 
with present and future fj, — e flavour experiments. Furthermore, an interesting prediction is 
the coincidence at the percent level of the conversion rate for Titanium and Lead experiments, 
in the full range of seesaw masses under 50 GeV. 

Considering our results and all constraints together, it follows that planned \x — > e con- 
version experiments will be fully relevant to detect or constrain sterile neutrino scenarios in 
an impressive mass range of 2 GeV - 1000 TeV. Note that this result is of general interest, as it 
may apply to any renormalizable BSM or dark matter theory which contains singlet fermions 
in its spectrum. 
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A Appendix: amplitudes 

As discussed in Section 3, the calculation of the [i — > e conversion rate requires the separation 
of the local and the long-range contributions. The latter is given by the dipole term in 
Eq. (3.1), while the local contribution is determined from Z and W mediated diagrams, as 
well as from the photon monopole interaction with the leptons. The amplitudes associated to 
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these three ingredients can be written as 15 



iM q , 

"/—monopole 



iM% 



2M, 



i a 



2M, 



w 



(A.l) 
(A.2) 
(A.3) 



In terms of the above form factors, F u and Fa as defined in Eq. (3.3) are given by 



(A.4) 
(A.5) 



In the case of light nuclei, the contribution of the dipolar term can also be casted in a four 
fermion interaction amplitude of the form 



ia w a_ „ 



W 



(A.6) 



In this way Eq. (3.5) is obtained: 



(A.7) 



The loop functions entering the computation of the form factors in Eqs. (A.4-A.6) are 



3 The term proportional to g M in Eq. (3.1) drops because of quark current conservation. 
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In terms of these functions, the form factors read 
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In the above, x 1)2 ,3 = a?n >2 , 3 = m ^, 2 , 3 / M W' x 4,...,3+fc = ^i,..., fe = m N li ... th / M W' x i = m g/ M w> 
V is the quark CKM matrix and U is the total (3 + k) x (3 + fc) neutrino mixing matrix 
defined in Eq. (2.3). The second equality in Eqs. (A.24)-(A.33) is obtained using the unitarity 
identity Yli^eiU^ = 0, the diagonalization relation ^j^Uei^/xlU^i = 0, the limiting values 
of the loop funtions in Eqs. (A.16)-(A.23) and the very good approximation for the present 
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analysis in which light neutrino masses are neglected compared to the heavy neutrino masses 
i x N h ... }k > £1/1,2,3)- 

These form factors present the following behaviour for low and high masses 
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where contributions involving 4 insertions of light-heavy mixing elements have been neglected; 
this is a good approximation in view of the various experimental bounds which hold on the 
mixings for the low energy regime, whereas for very high right-handed masses this approxi- 
mation relies on the perturbativity of the yukawa couplings. Dropping these terms with four 
insertions, the functions F q and F q can be written as 
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which are the expressions we actually use for all plots. 
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